Entangling strings of neutral atoms in ID atomic pipeline structures 
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We study a string of neutral atoms with nearest neighbor interaction in a ID beam splitter 
configuration, where the longitudinal motion is controlled by a moving optical lattice potential. 
The dynamics of the atoms crossing the beam splitter maps to a ID spin model with controllable 
time dependent parameters, which allows the creation of maximally entangled states of atoms by 
crossing a quantum phase transition. Furthermore, we show that this system realizes protected 
quantum memory, and we discuss the implementation of one- and two-qubit gates in this setup. 
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The recent development of optical and magnetic nii- 
crotraps allows the confinement of cold atoms in effective 
ID "pipeline" geometries, where the transverse quantum 
motion is frozen out . Variants of these trap designs 
promise the realization of beam splitters, and thus atomic 
interferometry "on a chip". Usually one envisions that 
atoms are injected one by one into these pipelines, where 
the source of cold atoms is provided by a Bose Einstein 
condensate. Instead we will study below collective beam 
splitter setups which allows the generation of entangled 
strings of atoms in ID trapping configurations with ap- 
plications in interferometry and quantum computing. 

To this end, we assume that the longitudinal motion of 
the atoms is controlled by storing atoms in a ID optical 
lattice potential generated by a standing light laser field. 
In the transverse direction the particles are confined by 
a a double well potential (see Fig. QJi) where we assume 
that the optical lattice stores exactly one atom per lattice 
site (i.e. one atom per double well). The preparation of 
such a Mott insulating state has been reported in a recent 
experiment, by loading of atoms from a Bose Einstein 
condensate via a superfluid - Mott insulator quantum 
phase transition . This setup by itself is an interesting 
extension of the standard "interferometry on a chip" , as 
it eliminates collisional shifts since atoms stored on dif- 
ferent lattice sites never collide. Furthermore, the atoms 
are supposed to be initially in the ground state which is 
a spatial superposition of the particles in the two trans- 
verse wells (region (I) of Fig. H^)- By moving the lattice 
we can drag the atomic chain "by hand" across the beam 
splitter while we increase the distance between the trans- 
verse wells adiabatically depending on the position of the 
atoms (i.e. we decrease the tunneling between the 
wells, see region (II) of Fig. ^). On the other hand, the 
use of optical lattices allows the engineering of coherent 
interactions between adjacent atoms (nearest neighbor 
interaction W in Fig. ^). This can be obtained either 
by cold collisions and moving optical lattices |^ j^, 
the remarkably strong dipole-dipole couplings of laser ex- 
cited Rydberg atoms Q, or by dipole-dipole coupfing of 
cold heteronuclear molecules 0- Together with appro- 
priate detection methods like fluorescence imaging these 




FIG. 1: External beam splitter: a) Atoms before (I) and after 
(II) the separation. The nearest neighbor interaction is de- 
noted by W and is the hopping matrix element between 
the two states \a) and |6) of the transverse trapping poten- 
tial. Internal beam splitter: b) atoms in two different internal 
states jo) and \b} enter the beam splitter. The internal states 
I a) and \b) are coupled by a Raman transition [see c)] with 
a Rabi frequency = fl. A laser excited Rydberg state |r) 
realizes the offsite interaction W with w denoting the width 
of the interaction zone. 



controllable interactions provide us with the tools to gen- 
erate entanglement of the ID chain of atoms. 

We will study the dynamics of the beam splitter setup 
indicated in Fig.^ In particular, we will investigate (i) 
how to generate a maximally entangled state of atoms, 
(ii) We will establish the formal equivalence of our model 
with well-studied models of spin chains. In particular, 
we will show that the system dynamics is a physical re- 
alization of a textbook model of a quantum phase tran- 
sition with completely controllable (time dependent) pa- 
rameters [^. Thus our setup provides an example of 
engineering a maximally entangled state from a prod- 
uct state via a quantum phase transition, (iii) Finally, 
the present setup implements the spin analogue 0| of 
Kitaev's protected quantum memory where qubits 
are represented by Majorana fermions, which provide a 
stable way to store quantum information due to an exci- 
tation gap 01 . Our setup allows to perform single and 
(collectively enhanced) two qubit operations. 
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We consider a ID chain of A'' atoms with modes \a) and 
1 5) stored in an optical lattice with a lattice constant A/2 
determined by the wave length A of the laser. The modes 
correspond either to two spatial modes in a double well 
structure, where the tunnelling provides a coupling (ex- 
ternal beam splitter in Fig. ^) , or to two internal atomic 
states connected via a Raman process (cf. Fig.n]D,c). We 
suppress hopping of the atoms between adjacent lattice 
sites by a sufficiently large potential barrier. This leads 
to an onsite interaction U ^ oo and we assume to have 
commensurate filling of one particle per lattice site. Fol- 
lowing U ^1 we derive a Hubbard Hamiltonian 

N-l 

H{t) = 2 ^ Wi{t) (aj^^ai+iajai + bj^^bi+iblbi'^ 
1=1 

N 

1=1 

Here describes coupling between \a) and \b) while the 
operators a;, 6/ are bosonic annihilation operators for 
particles in these modes at site I with [a;,6j] = [ai,6j] — 
0. A term J[ emerges from an additional state depen- 
dent superimposed trapping potential. We introduce the 
spin notation af — ajbi + aibj, — a\ai — bjbi and 
af = i{aibl — ajbi) which for ni = ajai + bjbi = 1 
are Pauli operators and rewrite the Hamiltonian as 

Hs{t) = EfjiWi{t)afaf,,-j:liJi{t) ■ <^i- Thus our 
setup is formally equivalent to an Ising chain of spins 
in a magnetic field J; = (Jf , Ji ,J[) |l3L.14]. 

Entanglement via quantum phase transition: Moving 
a string of atoms from left to right in the setup of 
Fig. ^i, or switching the lasers in Fig.QJjjC amounts to a 
time dependent change of the parameters from the large 
tunneling limit J^{t = 0) ^ \Wi\ to small tunnelling 
J^{t = T) ^ 0. In the following we assume that Jf'^ = 
except it is stated differently. In the homogeneous case 
(i.e. Jf = , Wi = W) the variation of amounts 
to crossing the critical point at = Vl^ of a quantum 
phase transition |^. Assuming that the atoms are ini- 
tially prepared in the product state | IT ' ' ' T)"^ with 
I T)r ~ + \b)i a superposition state of the two modes 
which is for = the (paramagnetic) ground state of 
Hs- Under adiabatic variation of parameters the system 
will remain in the ground state and evolve according to 
{W < 0) 

l+>= I !?•••?>" (ITT---T>^ + Ui---i>")/A/2 

^ (|0> + |1»/V2, (2) 

where the states | t)f = W)i and | |)f \b)i correspond 
to the atoms being in the upper or lower branch of the 
beam splitter of Fig. ^ The states |0) and |1) are the 
two degenerate (ferromagnetic) ground states of the Hs 
for = with all atoms in either one or the other arm 
of the beam splitter (see Fig.QJ. Thus the initial product 



state is transformed to a maximally entangled state via a 
quantum phase transition. The intuitive physical picture 
behind Q is as follows. Consider atoms moving across 
the beam splitter one by one. The first atom of the string 
will end up in the state | t)^ + | i)^) and attract the second 
atom. This leads to a state of the form | tt)^ + | [[ 
)^ . After the last atom has left the interaction zone the 
maximally entangled state |0) -I- |1) has been created. 

In the following we discuss the validity of the adia- 
batic approximation (Eq. ^) and thus the usefulness 
of this scheme by studying the scaling of the fidelity 
F — |(i/'id|'0(T'))P as a function of the length of the 
string TV and the time variation of Jf [t) and Wi (t) . Here 
F compares the state \il){T)) obtained from a time de- 
pendent integration of the Schrodinger equation with the 
ideal state li/iid) ^ |0) + |1)- This will be first done numer- 
ically, followed by analytical calculations and estimates. 

Before entering the time dependent case, we note that 
for the time independent case the Hamiltonian H has 
been studied extensively 0, 0| . For Jf = it can be 
fermionized and one obtains Hp = ^^^i^iftfu — 1/2) 
with the elementary excitation energies and fermionic 
annihilation (creation) operators (fl). The spectrum 
for the homogeneous case is shown in Fig. 12^.. For large 
iV the spectrum of the elementary excitations is charac- 
terized by a gap A = 2|W^ — for the energetically low 
lying quasi particles with the exception (arising from the 
free end boundary conditions) that the first excited state 
becomes degenerate with the ground or vacuum state 
(here, we do not take into account the second term in 
Hp, i.e. the vacuum state has zero energy) for \ W\ ^ 
(cf . Fig. |21i) . For — the two cat type ground states 
|0) -I- |1) and |0) — |1) correspond to the vacuum and the 
first excited state of the fermionized system, respectively. 

In Fig. we plot the numerically calculated operation 
time T required to perform (O with a fidelity of = 95% 
for linearly changing the homogeneous couplings J^{t) 
against N {W = const.). For iV > 20 we find a (poly- 
nomial) scaling of T for a given infidelity 1 — F ^ N"^ 
(cf. Fig. |2|3) in agreement with the analytical results be- 
low. By optimizing the time dependence of {t) we can 
speed up by the entanglement process significantly. 

A discussion of the spatially inhomogeneous situation 
where Jf and Wi vary as a function of I corresponding 
closer to the setup of Fig. ^ is given in Figs. |21;,d. For 
increasing time the string is moved across a zone of non 
vanishing Wi with a maximum and a width w. Si- 
multaneously, Jf (t) , is decreased from the initial value to 
Jf{T) « over a comparable "width" as w for all sites. 
The corresponding instantaneous time dependent energy 
levels are shown in Fig. [2t. Following the lowest energy 
curve in this diagram adiabatically from (1) to (2) cor- 
responds to |-|-) |0) -I- Fig. Eli shows the infidelity 
1 — for finite sweeping speed v against A'' for different 
widths w of the interaction zone. The infidelity 1 — F 
decreases rapidly with increasing w and scales exponen- 
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FIG. 2: Homogeneous setup Jf = , Wi = W: a) Elemen- 
tary excitations for = 25 against /W. The dashed 
horizontal line indicates the ground state energy, b) Upper 
bound F2 and lower bound Fi for the time T yielding a fidelity 
of -F = 95% as a function of A'^. The dashed line illustrates 
the A'^^ scaling predicted analytically. Inhomogeneous setup 
Ji , Wf. c) Instantaneous eigenenergies En for A'' = 6. d) 
Upper and lower bounds for the infidelity 1 — F against A' 
for constant sweeping speed v = O.OlATy and different inter- 
action zone widths w = O.IA (diamonds), w = 0.2A (closed 
circles) and w = 0.4A (open circles) and similar for . 



tially with for w < XN. For w > NX/2 the above 
scaling 1 — F ^ N"^ is restored. 

The numerical calculations behind Fig.|2]were based on 
a time independent Jordan- Wigner transformation of Hg, 
yielding a quadratic Hamiltonian in fermionic destruction 
(creation) operators {"I'D- By introducing Majorana 
operators C2u-i = {iv +7i)/2, = (71. - ll)l(2i) 
we obtain Ht — «c-^A(t)c, where the components of c 
are the Majorana operators and A(t) is a 27V x 2N real 
antisymmetric tridiagonal matrix. The linear Heisenberg 
equations of motion, c = A(t)c, are then solved numer- 
ically. We note that diagonalizing A in the time inde- 
pendent case yields Hp- For the fidelity F we use an 
approximate expression which can be derived as follows: 
The state |0) -I- |1) is the vacuum state of the fermion- 
ized system &t t ~ T. The completeness relation yields 
F{T) = l-E„^ol(n|^(r))P. Here |n) = |ni,...,n^) 
with Uy = 0, 1 the occupation numbers of the instan- 
taneous eigenstates of Ht corresponding to an energy 
ty. The sum in this expression can be reordered, and 
we obtain F{T) = 1 — where P{m) = 

X^n I ('^IV'(T')) P'^m.j; m is the probability of having m ele- 
mentary excitations in the system at time T . By solving 
the above equation for c we can in principle calculate 



the quantities Ai = flfm)^) = P{m)mK 

The fidelity F is then given by the solution of a system 
of N linear equations. An approximate fidelity Fi can be 
obtained by neglecting the probabilities P{k) with k > I. 
We restrict ourselves to / = 1, 2 and find Fi — 1 — Ai and 
F2 = 1 - (3^1 - ^2)72. The exact fidelity is bounded 
by these quantities: Fi < F < F2. Compared to a cal- 
culation in the spin picture which requires the solution 
of ^ 2^ equations the calculation of Ai and A2 can be 
done by solving a system of ~ A'^^ differential equations. 

Let us turn to the more technical point of analyti- 
cally estimating the scaling of the fidelity F when the 
phase transition point is crossed by linearly changing 
— Qt + W with Q =const. First we note that there 
are no transitions between the ground and the first ex- 
cited state since they have opposite parity. Close to the 
phase transition point the energy gap to the remaining 
excitations A « and therefore at the time t = — the 
evolution of the system ceases to be adiabatic and excita- 
tions start to be populated. The adiabaticity is restored 
again at the time t t^, when the gap A becomes suf- 
ficiently large to prevent further excitations. Then, the 
relaxation of the new phase occurs separately within dif- 
ferent domains, whose sizes are given by the value lo{t*) 
of the correlation length at the time t*. Close to the 
phase transition Iq ^ A~^/^ and therefore the domain 
size scales like /o(^*) ^ O^^^^. The quench through the 
phase transition point can only be adiabatic if the char- 
acteristic size of the domain formed exceeds the size of 
the system L ^ N and therefore ^o(^*) ?J L, which gives 
the scaling condition 6 < W'^/N'^, or WT ~ iV^. 

Quantum computing model with protected quantum 
memory: In the case < the ferromagnetic super- 
position state is very sensitive to homogeneous distor- 
tions of the form J[ = which induce a relative phase 
shift exp(j2iV /J" dtJ^t)) scaling with N ^ between the 
two states |0) and |1) after a time r. Therefore, in the 
external beam splitter setup where these two states are 
spatially separated they can be viewed as two arms of a 
Heisenberg limited interferometer collectively enhanced 
by a factor N. On the other hand, in the antiferromag- 
netic case, i.e. for a repulsive interaction W > 0, the two 
degenerate ground states at = 



|o) = uT•••iT)^ |l) = lTi•••n)^ 



(3) 



are closely related to unpaired Majorana fermions which 
have been considered as candidates for storing quantum 
information 0, 0|. These states are expected to be in- 
sensitive against perturbations since they are separated 
by a gap of order W from the other states of the sys- 
tem and are only connected via iV-th order perturbation 
theory for homogeneous couplings 3i — 3. This yields 
stability against spin flip errors exponentially increasing 
with the number of particles in the chain N and is also 
reflected by the scaling of the energy of the first excited 
state ei - (J^/l^)^ for < 0. Furthermore, 
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a) ,,.^j/^^A^,_,,b) ^ .yy^,. 
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FIG. 3: Collectively enhanced interactions between two 
strings of atoms 1 and 2. a) Antiferromagnetic setup: N/2 
particles of each chain interact with strength W' only if they 
are in states |0)i|l)2 or |l)i|0)2 yielding a phase gate between 
the two qubits implemented by those chains, b) Ferromag- 
netic setup: Entanglement creation between two chains of 
atoms via interactions W' in the state |l)i|0)2. Impelemen- 
tations with optical lattices or atom chips, for instance, offer 
the scalability of the scheme. 
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FIG. 4: Illustration of the Hadamard gate for = 8 by 
adiabatically changing and (unprotecting the quantum 
memory) . We follow the lowest two eigenstates (with energies 
given by the solid curves) transforming as |0) + |1) — > |0) and 
|0) — |1) — > |1) (up to a dynamical phase) when changing J^, 
in three steps (1), (2), (3) (followed by turning off in 
step (4)) as described in the text. Note that if the condition 
J' < W/{N—1) is not fulfilled we get unwanted crossings and 
the first excited state after step (3) will not be of the form 
The dashed curve shows the third eigenenergy of the system 
and the inset the path in the — plane. 



if we assume that N is even the states |0) and |1) are 
completely insensitive to global fluctuations of J since 
EiO-f|0> = EiO-ill) = 0- Then the two states |0) and 
|1) constitute a decoherence free subspace 0,^^ for ho- 
mogeneous perturbations and can thus be used as qubits 
which store quantum information reUably. 

We will now discuss how to implement single and (col- 
lectively enhanced) two qubit gates and show that our 
model realizes a quantum computer with protected mem- 
ory. The idea behind the two qubit phase gate is sum- 
marized in Fig. 13 Selectively overlapping the wave func- 
tions of different two qubit states for a time T2 colli- 
sional interactions of strength W' between the atoms 
lead to an entanglement phase <p2 — NW't2/2 Q cor- 
responding to a phase gate with a truth table |ei)|e2) ^ 
exp(i02((ei + e2)mod2))|ei)|e2), with ei^2 = {0,1}. Sin- 
gle qubit gates correspond to a general unitary rotation 
of |0) and |1) (Eq. Q) which can be decomposed in 
Hadamard gates |0), |1) |0) ± |1) and the creation of a 



relative phase |ei) — > exp(iei(/)i)|ei). The phase (pi can 
be implemented by turning on a trap potential creating 
a staggered offset of the form — J^(— 1)' for a time 
Ti = (pi/2NJ^. The idea behind the Hadamard gate is 
as follows: At J"^ = = the states |0), |1) represent a 
degenerate eigenspace of Hs- Turning on the field up 
to > W, thus unprotecting the qubit and switching it 
off when ^ will under appropriate conditions induce 
a rotation in this space. A specific example is illustrated 
in Fig. 01 (1) at = we adiabatically switch on 
until > is reached unprotecting the qubit, then (2) 
we increase J^, (3) we return adiabatically to = 0, 
and, finally, (4) switch off J^. 

We have shown how to generate maximally entangled 
states of strings of atoms in ID pipeline configurations. 
An extension of this setup allows implementations of a 
quantum computing model with protected qubits. 
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